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PREFACE 

The  research  reported  herein  was  performed  by  the  Arnold  Engineering 
Development  Center  (AEDC) , Air  Force  Systems  Command  (AFSC) , under 
Program  Element  6S807F.  Work  and  analysis  for  this  research  was  done  by 
personnel  of  ARO,  Inc. , AEDC  Division  (a  Sverdrup  Corporation  Company) , 
operating  contractor  of  AEDC,  AFSC,  Arnold  Air  Force  Station,  Tennessee. 
The  work  was  conducted  under  ARO  Project  No.  P32A-C1A.  The  author  of 
this  report  was  M.  C.  Altstatt,  ARO,  Inc.  The  manuscript  (ARO  Control 
No.  ARO-FWT-TR-77-22)  was  submitted  for  publication  on  March  28,  1977. 
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1.0  INTRODUCTION 

In  transonic  wind  tunnel  testing  of  airfoils,  the  Reynolds  number 
scaling  presents  a difficult  problem  because  of  the  complex  interaction 
of  the  shock  wave  with  the  boundary-  layer.  Generally,  flow  separation 
is  induced  at  or  near  the  shock  and  this  separation,  in  turn,  affects 
the  location  of  the  shock  wave.  In  order  to  predict  the  aerodynamic 
coefficients  of  a wing,  especially  the  moment  coefficient,  an  accurate 
prediction  of  the  shock  location  is  necessary.  The  difficulty  in  pre- 
dicting the  change  in  this  location  with  change  in  Reynolds  number 
reflects  an  inadequate  understanding  of  the  physical  phenomena  and 
therefore  in  the  equations  used  to  describe  these  phenomena. 

To  provide  the  basic  data  required  to  evaluate  certain  computa- 
tional methods,  two  wind  tunnel  tests  were  conducted.  The  data  obtained 
included  pitot  pressure  profiles  in  the  boundary  layer  and  laser  velocim- 
eter  (LV)  measurements  in  the  boundary  layer  and  free  stream.  A com- 
puter program  has  been  developed  to  reduce  the  LV  data,  with  emphasis  on 
detecting  and  minimizing  sources  of  measurement  biasing. 

Computational  results  for  a turbulent  compressible  boundary  layer, 
obtained  using  several  different  Reynolds  stress  models,  are  compared 
with  the  experimental  results.  The  calculations  are  not  accurate  where 
pressure  gradients  are  large  or  near  flow  separation.  However,  the 
models  used  for  Reynolds  stress  and  boundary- layer  velocity  produce 
results  which  do  compare  well  with  experimental  data,  if  based  on  measured 
boundary- layer  thickness  and  skin  friction. 
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2.0  EXPERIMENTAL  PROGRAM 

Two  tests  were  conducted  in  the  Aerodynamic  Wind  Tunnel  (IT).  A 
sketch  of  the  test  section  arrangement  is  shown  in  Fig.  1.  The  model  is 
a circular  arc  bump,  modified  by  a cosine  fairing  at  each  end,  mounted 
on  the  tunnel  bottom  wall.  Testing  was  conducted  at  a nominal  Mach 
number  of  0.8  and  a unit  Reynolds  number  of  approximately  five  million 
per  foot.  The  pressure  distribution  and  velocity,  as  measured  by  sur- 
face static  ports,  are  shown  in  Fig.  2a.  From  Fig.  2a  and  the  shadow- 
graph, Fig.  2b,  it  can  be  seen  that  this  is  a transonic  flow  of  the 
type  desired,  with  a shock-wave/boundary-layer  interaction  and  flow 
separation.  The  supercritical  region  extends  from  approximately  station 
18  to  21 , where  the  shock  occurs.  The  separation  point  is  at  about 
station  21.75. 

2.1  TESTING 

For  the  first  test,  22  pitot  pressure  profiles  were  obtained  in  the 
boundary  layer.  The  probes  were  projected  through  the  tunnel  bottom 
wall  and  positioned  by  an  electric  drive.  High-speed  shadowgraph  motion 
pictures  were  taken  to  determine  the  flow  stability.  At  400  frames  per 
second  the  boundary  layer,  shock  position,  and  separation  point  appeared 
to  be  stationary  on  the  surface,  but  at  a height  of  above  2 in.  the 
shock  was  oscillating  slightly. 

In  the  second  test,  a laser  velocimeter  system,  described  in  Ref. 

1,  was  used  to  measure  flow  velocity  profiles  at  19  boundary-layer 
profiles  and  four  streamwise  paths  through  the  shock.  The  program 
developed  for  the  LV  data  reduction  is  described  in  the  Appendix.  A 
streamwise  coordinate  system  is  used  for  all  the  LV  boundary- layer 
results  presented. 
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2.2  TEST  RESULTS 

For  direct  comparison  of  pitot  and  LV  velocity  measurements,  several 
boundary-layer  profiles  taken  with  the  pitot  probe  in  the  first  test 
were  repeated  with  the  LV  system  in  the  second  test.  Figure  3a  shows  a 
typical  comparison  for  data  taken  near  the  leading  edge  of  the  bump. 

There  is  good  agreement  for  all  profiles  taken  at  corresponding  locations 
in  the  attached  flow. 

In  the  separated  region,  as  illustrated  in  Fig.  3b,  the  LV  and 
pitot  probe  results  do  not  agree.  The  LV  data  are  preferred  in  the 
separated  region  because  the  pitot  probe  is  affected  by  variations  in 
flow  angle  and  possible  variations  in  local  static  pressure  and  free- 
stream  stagnation  pressure.  Sources  which  can  affect  the  interpretation 
of  LV  measurements,  such  as  particle  lag,  velocity  biasing,  and  time 
variation,  are  discussed  in  Refs.  2 and  3 and  in  the  Appendix. 

A problem  in  consistently  Interpreting  the  LV  data  arises  in  deter- 
mining the  location  of  the  edge  of  the  boundary  layer.  The  pitot  data 
give  consistent  criteria  for  the  edge,  as  long  as  the  free-stream  stag- 
nation pressure  remains  constant.  With  velocity  data  only,  however,  as 
canbe  seen  from  Fig.  3c,  the  boundary-layer  edge  is  less  clearly  defined. 
This  ambiguity  can  affect  the  computation  of  the  boundary- layer  properties. 

Starting  above  the  bump  centerline,  LV  velocity  surveys  were  con- 
ducted downstream  parallel  to  the  tunnel  bottom  wall.  Figure  4 shows 
velocity  and  stream  angle  for  four  profiles  at  1.25,  1.5,  1.75,  and  2.0 
in.  above  the  bottom  wall.  These  profiles  indicate  a variation  from  an 
almost  normal  shock  at  y - 2.0  in.  to  a family  of  oblique  shocks  at  the 
surface. 

The  LV  data  presented  in  this  report  were  processed  by  a data 
reduction  program  which  is  described  in  the  Appendix.  Information, 
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such  as  mode,  standard  deviation,  skewness,  kurtosis,  and  covariance, 
is  extracted  from  the  LV  velocity  measurements.  As  defined  in  Ref.  4, 
the  standard  deviation  and  covariance  are  directly  related  to  the  tur- 
bulent intensity  and  Reynolds  stress.  The  mode,  skewness,  and  kurtosis 
are  assessed  for  evidence  of  biasing  in  the  measurements. 

Mean  and  mode  velocity,  turbulent  intensity,  and  Reynolds  stress 
profiles  in  the  separated  flow  at  station  23.88  are  shown  in  Fig.  5. 
Because  of  the  flow  fluctuation  in  this  region  and  the  length  of  time 
required  to  obtain  these  measurements,  the  accuracy  may  be  poor  near 
the  surface. 


3.0  THEORETICAL  ANALYSIS 

Computations  were  performed  for  the  boundary  layer  measured,  using 
the  experimental  pressure  distribution  as  an  input.  The  solution  was 
compared  with  the  test  data  to  determine  the  accuracy  of  the  theoretical 
solutions  and  evaluate  available  turbulence  models. 

3.1  BOUNDARY-LAYER  CALCULATION  METHOD 

The  calculation  method  used  for  a turbulent,  compressible  boundary 
layer  is  a variation  of  the  method  of  Nash  and  Hicks  (Ref.  5)  as  modified 
by  Kuhn  and  Nielsen  (Ref.  6).  The  formulation  and  computer  code  used 
are  reported  in  Ref.  7. 

A summary  of  the  equations  is  given  below. 

Continuity 


/-  (P“)  + (Pv)  = 0 

OX  OY 


(D 
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Momentum  Energy 


Energy 


pu 


du  du  dP  d ( a du\ 

r,*PWTy-h  + iyl(fi^) 


pa 


dS  dS  d ( a ds\ 


(2) 


(3) 


It  is  assumed  that  the  Prandtl  number  is  unity  for  both  laminar  and 
turbulent  flow.  In  these  equations,  (3  is  defined  by 


/S  = 1 + tfv 


(4) 


where  t is  the  eddy  viscosity, 


pu  v 


(5) 


u'  and  v'  are  the  turbulent  velocity  fluctuations,  and  the  bar  denotes  a 
time  average.  The  enthalpy  parameter,  S,  is  given  by 


S = T/Tl - 1 

with  T as  the  total  temperature. 


(6) 


Application  of  the  Stewartson  transformations  to  these  equations 
results  in  the  set: 


dll  W 
dX  dY 


Momentum 


dr  dll  P dU 

U— +V -U  — - 

ax  av  p ” d\ 


(7) 


(8) 
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The  terms  U,  V,  X,  and  Y are  transformed,  incompressible  velocities, 
and  coordinates,  and  the  subscript  e refers  to  conditions  at  the  edge  of 
the  boundary  layer.  An  adiabatic  wall  is  assumed  and  Crocco’s  relation 
with  a recovery  factor  of  0.89  is  used  as  an  empirical  solution  to  the 
energy  equation. 


— = 1.0  + 0.178  M2 
P 


(9) 


Substituting  Eq.  (7)  into  Eq.  (8)  and  integrating  across  the  boundary 
layer  results  in 

if[UU,UY/IUI(d9-^UeUaI-I,^UY>Y]  YndV  = 0 <10> 

where  n = 0,  1 . 


A wall-wake  model  is  used  for  the  boundary-layer  velocity  profile, 

U = ur  [2.5  £«  (1  + y+)  + 5.1  -{3.39y+  + 5.1)(!  “°-375f+]  + ^^l  - cos  (11) 

where  U^.  is  a friction  velocity, 

V’'J  .02) 


Ug  is  a wake  velocity, 

up  = lJe  — Uy.  [2.5  (1  + 5*)  + 5.1] 


03) 


and 


y+  = ur  Y/ v 


(14) 


This  set  of  equations  can  be  closed  with  a model  for  eddy  viscosity. 
The  inner  and  outer  eddy  viscosity  models  given  by  Alber  (Ref.  8)  are 
used : 


Inner  wall  model 


/3,  = 1 .0  + 0.0533  teu+  - [1.0  + U+  T !,  IT2!  \ 


(15) 
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where 

U+  = 0.41  [}/uT  (16> 

Outer  wake  model 

US* 

0O  = 1.0  + 0.0168  [1.0  + 5.5  (Y/5)6]-1  (17) 

v 

A set  of  differential  equations  was  therefore  formed,  in  terms  of  u^ 
and  6,  and  the  Runge-Kutta  method  of  advancing  in  the  streamwise  direc- 
tion was  used. 


3.2  COMPARISON  WITH  EXPERIMENT 

Figure  6 Illustrates  measured  and  calculated  displacement  thickness 
on  the  model.  The  agreement,  as  shown  in  Fig.  6a,  is  good  to  about 
tunnel  station  16,  where  the  pressure  gradient  is  becoming  large.  Past 
that  point,  the  computed  displacement  thickness  deviates  rapidly  from 
the  experimental  values.  Figure  6b  shows  a comparison  of  analytic  and 
measured  displacement  thickness  downstream  in  the  region  of  the  shock 
and  separation.  Since  the  solution  shown  in  Fig.  6a  became  unstable 
past  station  16,  this  solution  was  restarted  at  station  18.  As  expected, 
the  boundary-layer  solution  is  not  accurate  in  this  region  of  strong 
pressure  gradient  and  flow  reversal. 

3.3  REYNOLDS  STRESS  MODELS 

In  order  to  determine  the  sensitivity  of  the  solutions  to  the 
modeling  of  the  Reynolds  stress  term  in  Eq.  (10),  several  expressions 
for  the  eddy  viscosity  were  used.  These  forms  are 


08) 
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Michael-Quemard-Durant  (Ref.  9) 


f = F2  £2^’ 
dY 


F = 1.0  - e 


V' 


0.0855 


£ = 0.085  8 Unh 


Lees  (modified)  (Ret.  8) 


(19) 

(20) 
(21) 


(22) 


A laminar  case  (e  = 0)  was  also  run.  In  the  cases  for  which  no  inner 
model  is  given,  Eq.  (15)  is  used. 

A comparison  of  the  resulting  Reynolds  stresses  (0)  is  given  in 
Fig.  7.  Although  these  models  show  substantial  differences  in  0 dis- 
tribution, the  effect  is  small  on  the  computation  of  displacement  thick- 
ness, for  the  region  near  the  bump  LE,  as  shown  in  Fig.  8.  Also,  down- 
stream of  station  16  the  solution  becomes  unstable  regardless  of  the 
eddy  viscosity  model.  The  local  variations  across  the  boundary  layer 
are  smoothed  by  the  Integration  process. 

1 

3.4  VERIFICATION  OF  ANALYTIC  MODELS 


Results  from  the  models  used  for  the  boundary-layer  velocity  profiles 
[Eq.  (11)]  and  Reynolds  stress  [Eqs.  (15)  and  (17)]  can  be  compared 
directly  with  experimental  measurements.  Figures  9a  and  b illustrate 
such  a comparison  at  tunnel  station  14.5.  Even  though  the  agreement 
between  the  calculated  and  experimental  values  shown  is  very  good,  it  is 
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at  about  this  point  where  the  downstream  advance  by  the  boundary- layer 
scheme  begins  to  deviate  sharply  from  the  measured  boundary  layer.  This 
deviation  may  exist  because  the  momentum  balance  for  the  streamwlse 
advance  is  inappropriate  in  the  large  pressure  gradient. 

As  noted  above,  the  boundary-layer  9cheme  is  not  valid  in  the 
shock— separation  region  downstream.  However,  at  locations  where  experi- 
mental data  are  available,  Eqs.  (11)  through  (17)  can  be  used  with  the 
measured,  integrated  boundary- layer  properties  as  input,  and  velocity 
profiles. and  Reynolds  stress  computed,  as  in  starting  solutions. 

Figures  9c  and  d show  a comparison  of  these  computations  with  the 
measured  velocity  and  Reynolds  stress.  It  can  be  seen  that  even  though 
the  calculation  is  not  valid  in  the  downstream  direction,  the  models 
UBed  in  the  equations  are  consistent  with  experimental  data  at  a given 
station.  Therefore,  the  boundary-layer  equations  utilized  are  suspect 
in  this  region. 


4.0  CONCLUSIONS 

An  investigation  of  a transonic  flow  containing  a shock-wave/boundary- 
layer  interaction  of  the  type  occurring  on  airfoils  was  conducted,  using 
both  experimental  and  theoretical  methods.  Tests  were  conducted  in  the 
Aerodynamic  Wind  Tunnel  (IT).  Data  Included  surface  static  pressures, 
high-speed  shadowgraph  motion  pictures,  and  pitot  pressure  profiles  in 
the  boundary  layer.  A laser  velocimeter  system  was  used  to  obtain 
velocity  measurements  in  the  boundary  layer  and  in  the  free  stream.  A 
data  reduction  program  was  developed  to  extract  mode  and  mean  velocity, 
standard  deviation,  skewness,  kurtosis,  and  Reynolds  stress  from  the 
velocity  measurements. 

A program  for  a turbulent,  compressible  boundary  layer  was  used  to 
compute  the  boundary- layer  properties  for  the  flow  field  investigated. 

A comparison  of  the  analytic  and  experimental  displacement  thickness 
shows  good  agreement  in  attached  flows  except  in  regions  with  large 
pressure  gradients  and  near  the  separation  point. 
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An  investigation  was  initiated  to  determine  the  accuracy  of  some  of 
the  turbulence  models  used.  Several  eddy  viscosities  were  used  in  the 
computation  of  Reynolds  stress  for  the  boundary- layer  equations.  Even 
though  the  various  Reynolds  stresses  calculated  were  substantially 
different,  the  resulting  differences  in  displacement  thickness  were 
relatively  small.  The  integral  boundary-layer  scheme  smooths  out  the 
local  variations  in  eddy  viscosity. 

i 

A comparison  of  the  computed  and  measured  Reynolds  stresses  and 
0 

boundary-layer  velocity  profiles  at  selected  points  shows  good  agree- 
ment, even  where  the  computed  boundary- layer  results  are  deviating  from 
the  measured  values.  The  set  of  boundary- layer  equations  used  is  appar- 
ently inadequate  for  the  regions  of  strong  pressure  gradient. 

In  the  separated  region  to  the  rear  of  the  shock,  good  agreement 
can  also  be  obtained  between  the  computed  and  measured  Reynolds  stresses 
and  velocities  if  experimental  data  are  used  as  input  in  the  calculation. 
Even  though  the  method  of  solution  for  the  boundary  layer  is  inaccurate, 
the  turbulence  models  are  self-consistent  at  any  given  station. 
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a.  Velocities  upstream  of  the  model,  station  9.35 
Figure  3.  Comparison  of  laser  velocimeter  and  pitot  tube 
measurements  in  the  boundary  layer. 
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b.  Velocities  in  separated  flow,  stetion  23.88 
Figure  3.  Continued.  ' 
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a.  Velocity  profile 


b.  Turbulent  intensity 


c.  Reynolds  stress 


Figure  5.  Boundary-layer  properties. 
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TUNNEL  STATION,  inches 

b.  Shock  wave/boundary-layer  interaction  and 
separation  regions 

Figure  6.  Concluded. 
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a.  Velocity  profile  at  station  14.5 
Figure  9.  Comparison  of  analytic  results  with  experimental 
measurements. 
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b.  Reynolds  stress  at  station  14.5 
Figure  9.  Continued. 
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c.  Velocity  profile  at  station  23.88 
Figure  9.  Continued. 
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d.  Reynolds  stress  at  station  23.88 
Figure  9.  Concluded. 
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The  laser  velocimeter  (LV)  system  records  two  components  of  velocity 
for  a particle  passing  through  the  beam  intersection  volume,  which  i9 
approximately  0.3  mm  in  diameter  by  9 mm  in  length.  More  details  of  the 
LV  and  data  acquisition  system  are  reported  in  Ref.  10.  At  a point  in 
the  flow,  approximately  1,000  velocity  pairs  are  taken.  The  time  required 
to  obtain  the  data  can  vary  from  leas  than  a second  to  several  minutes, 
depending  on  the  number  of  particles  per  unit  time  passing  .through  the 
flow  volume.  Figure  A-1  illustrates  a data  set  which  is  low  in  the 
boundary  layer  near  the  separation  point. 

DATA  REDUCTION  EQUATIONS 

The  flow  characteristics  computed  from  the  LV  data  are  mean  and 
mode  velocity,  turbulent  intensity,  skewness,  kurtosis,  and  Reynolds 
stress.  The  following  expressions  were  used: 

Mean  velocity 


Standard  deviation 


Skewness 


Kurtosis 


N 


■ - i £ 

i=i 


(A-1) 


N i=l 


(A-2) 


(uj  - u)3 
N<7  3 1 


(A-3) 


Nct4  1 
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Convarience 

, « 

u'v'  = Vs  (u.  - u)(v.  - v)  (A-5) 

N . 

]=  1 

where  and  are  Individual  velocity  samples  and  N is  the  number  of 
samples . 

To  determine  the  mode  velocity,  a histogram  is  formed,  as  in  Fig. 
A-2.  Here  the  number  of  velocity  samples  falling  within  a given  velocity 
range,  usually  10  ft/sec,  is  plotted  against  velocity.  A least-squares 
polynominal  curve  fit  (the  curve  shown  on  the  plot)  is  obtained,  and  the 
maximum  of  the  curve  is  taken  as  the  mode  value. 

DATA  INTERPRETATION 

Although  the  data  reduction  equations  are  simple,  interpretation  of 
the  results  presents  problems.  Some  of  the  factors  producing  uncer- 
tainty in  the  results  are  listed, below. 

Tracking:  In  the  unseeded  flow  used,  the  particle  size  and  composition 

distribution  is  not  known,  so  there  can  be  an  unknown  difference 
between  particle  velocity  and  fluid  velocity  for  the  measurements. 

Velocity  Bias:  For  some  combinations  of  velocity  and  particle  number 

density,  as  shown  in  Ref.  3,  the  probability  of  obtaining  high- 
velocity  data  is  greater,  since  more  high-velocity  particles  pass 
through  the  probe  volume. 

Flow  Fluctuation:  Because  of  the  time  required  to  obtain  1,000  velocity 

sets,  if  the  flow  is  fluctuating  or  tunnel  conditions  are  varying. 
Improper  readings  are  obtained. 
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Position  Drift;  Especially  for  low  boundary- layer  points,  where  posir 
tion  is  critical  and  record  times  can  be  long,  temperature  change 
and  vibration  can  produce  error. 

Processor  Bias:  The  LV  signal  processor  is  biased  toward  certain  digits 

in  the  last  decimal  place,  which  introduces  an  error  if  the  stand- 
ard deviation  is  small. 

The  uncertainty  Introduced  by  these  factors  can  vary  from  less  than 

i ' ‘ i i ■ 

2 percent  in  a free-stream  flow  to  the  order  of  20  percent  for  a boundary- 
layer  point  in  a fluctuating  flow.  Some  of  the  options  built  into  the 
data  reduction  program  are  designed  to  indicate  or  minimize  these  errors. 

PROCESSING  OPTIONS 

The  data  processing  options  in  the  program  fall  into  three  cate- 
gories: selection  of  valid  data,  method  of  processing,  and  presentation 
of  results. 

Given  a set  of  1,000  two-component  velocity  readings,  some  of  these 
will  have  one  or  both  measurements  misread.  These  bad  readings  are 
removed  from  the  data  set.  Then,  the  data  set  can  be  further  reduced  by 
the  deletion  of  velocity  records  which  differ  from  the  mean  by  more  than 
some  value,  usually  three  standard  deviations.  Also,  the  data  set  can 
be  reduced  to  simultaneous  components  or  all  velocities  of  each  compo- 
nent can  be  processed. 

The  output  can  be  presented  in  tunnel  or  streamwise  coordinates  or 
rotated  to  any  desired  angle  from  either  system.  Print  options  vary 
from  a one-line  summary  at  each  point  to  a full  print  with  raw  data  and 
histograms. 

Variations  in  processing,  checking  for  extreme  values  of  the  velocity 
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moments,  comparing  mean  and  mode  velocities,  and  other  indicators  from 
the  summary  outputs  determine  which  data  points  require  additional 
study.  In  some  cases,  errors  can  be  identified  and  corrected,  or  the 
data  point  is  discarded. 
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Figure  A-1.  Data  set  near  the  separation  point. 
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Figure  A-2.  Histogram  showing  mode  velocity. 
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NOMENCLATURE 

Mixing  length 
Number  of  samples 
Pressure 
Reynolds  number 
Enthalpy  parameter 
Temperature 

Transformed,  incompressible  velocities 
Mean1  velocities 
Instantaneous  velocities 
Turbulent  velocity  fluctuations 
0.041  U/ut 

Wake  velocity  parameter 
Friction  velocity  parameter 
Transformed,  incompressible  coordinates 

Physical  coordinates;  Y is  measured  from  surface  for  boundary- 
layer  flow,  above  tunnel  floor  otherwise;  x is  streamwise  for 
boundary  layers. 
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Coordinate  at  edge  of  sublayer 

y+  uTy/v 

6 1 + e/v 

5 Boundary-layer  thickness 

6*  Displacement  thickness 

e Eddy  viscosity 

6 Momentum  thickness 

v Kinematic  viscosity 

V Viscosity 

p Density 

a Standard  deviation 

SUBSCRIPTS 

AW  Adiabatic  wall 

e Boundary-layer  edge 

i Inner 

o Outer 

t Stagnation 

x,y  Partial  derivative 

6 Referred  to  boundary-layer  thickness 
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